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Generating Functions For Kernels of Digraphs 
(Enumeration & Asymptotics for Nim Games) 

Cyril Banderier, Jean-Marie Le Bars, and Vlady Ravelomanana 

Abstract. In this article, we study directed graphs (digraphs) with a coloring constraint due 
to Von Neumann and related to Nim-type games. This is equivalent to the notion of kernels of 
digraphs, which appears in numerous fields of research such as game theory, complexity theory, 
artificial intelligence (default logic, argumentation in multi-agent systems), 0-1 laws in monadic 
second order logic, combinatorics (perfect graphs)... Kernels of digraphs lead to numerous diffi- 
cult questions (in the sense of NP-complotencss, #:P-completeness). However, we show here that 
it is possible to use a generating function approach to get new informations: we use technique 
of symbolic and analytic combinatorics (generating functions and their singularities) in order to 
get exact and asymptotic results, e.g. for the existence of a kernel in a circuit or in a unicircuit 
digraph. This is a first step toward a generatingfunctionology treatment of kernels, while using, 
e.g., an approach "a la Wright". Our method could be applied to more general "local coloring 
constraints" in decomposable combinatorial structures. 

Resume. Nous etudions dans cet article les graphes diriges (digraphes) avec une contrainte de 
coloriage introduite par Von Neumann et reliee aux jeux de type Nim. EUe equivaut a la no- 
tion de noyaux de digraphes, qui apparait dans de nombreux domaines, tels la theorie des jeux, 
la theorie de la complexite, I'intelligence artificielle (logique des defauts, argumentation dans 
les systemes multi-agents), les lois 0-1 en logique monadique du second ordre, la combinatoire 
(graphes parfaits)... Les noyaux des digraphes posent de nombreuses questions difficiles (au sens 
de la NP-completude ou de la #P-completude) . Cependant, nous montrons ici qu'il est possible 
de recourir aux series generatrices afin d'obtenir de nouvelles informations : nous utilisons les 
techniques de la combinatoire symbolique et analytique (etude des singularites d'une serie) afin 
d'obtenir des resultats exacts ou asymptotiques, par exemple pour I'existence d'un noyau dans 
un digraphe unicircuit. II s'agit la de la premiere etape vers une serie generatrilogie des noyaux. 
Notre methodc peut etre appliquee plus generalement a des "contraintes locales" de coloriage 
dans des structures combinatoires decomposables. 



1. Introduction 

Let V and E be the set of vertices and directed edges (also called arcs) of a directed graph D 
without loops or multiarcs (we call such graphs digraphs hereafter) . A kernel of D is a nonempty 
subset K of V, such that for any a,b £ K, the edge (a, b) does not belong to E, and for any vertex 
outside the kernel (a ^ K), there is a vertex in the kernel (6 G K), such that the edge (a, b) belongs 
to E. In other words, X is a nonempty independent and dominating set of vertices in D [2]. Not 
every digraph has a kernel and if a digraph has a kernel, this kernel is not necessarily unique. The 
notion of kernel allows elegant interpretations in various contexts, since it is related to other well- 
known concepts from graph theory and complexity theory. In game theory the existence of a kernel 
corresponds to a winning strategy in two players for famous Nim-type games (cf. [3, 16, 17, 31]). 

Imagine that two players A and B play the following game on D in which they move a token 
each in turn: A starts the game by choosing an initial vertex vq G V and then makes a move to 
a vertex vi. A move consists in taking the token from the present position Vi and placing it on a 
child of Vi, i.e. a vertex Vi+i such that (vi,Vi+i) ^ E. B makes a move from vi to V2 and gives 
the hand to A, which has now to play from U2, and so on. The first player unable to move loses 
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the game. One of the two players has a winning strategy (as this game is finite in a digraph D 
without circuit, for circuits one extends the rules by saying that the game is lost for the player 
who replays a position previously reached). Von Neumann and Morgcnstcrn [31] proved that any 
directed acyclic graph has a unique kernel, which is the set of winning positions for A {A always 
forces B to play outside the kernel, until B cannot play anymore). Richardson [27] proved later 
that every digraph without odd circuit has a kernel [7, 29]. Berge wrote a chapter on kernels 
in [2] . Furthermore, there is a strong connection between perfect graphs and kernels (see the Berge 
and Duchet survey [1]). Some natural variants of this property are studied in various logic for 
Intelligence Artificial, some of them are definable in default logic [8] and used for argumentation 
in multi-agents systems, kernels appear there as sets of coherent arguments [6, 12]. 

Fernandez de la Vega [13] and Tomescu [30] proved independently that dense random digraphs 
with n vertices and m = <d{'n?) edges, have asymptotically almost surely a kernel. In addition, 
they get the few possible sizes of a kernel and a precise estimation of the numbers of kernels. 

Few years ago a new interest for these studies arises by their applications in finite model 
theory. Indeed variants of kernel are the best properties to provide counterexamples of 0-1 laws in 
fragments of monadic second-order logic [21, 22]. Goranko and Kapron showed in [19] that such 
a variant is expressible in modal logic over almost all finite frames for frame satisfiability; recently 
Le Bars proved in [23] that the 0-1 law fails for this logic. 

The existence of a kernel in a digraph has been shown NP-complete, even if one restricts this 
question to planar graphs with in- and out-degree < 2 and degree < 3 [9, 11, 15]. It is somehow 
related to finding a maximum clique in graphs [4, 21], which is known to be difiicult for random 
dense graphs. 

In this article, we use some generating function techniques to give some new results on Nim-type 
games played on directed graphs (or, equivalcntly, some new informations on kernel of digraphs). 
More precisely, we deal with a family of planar digraphs with at most one circuit or one cycle and we 
give enumerative (Theorems 4.1, 4.2, 4.3, 4.4 in Section 4) and asymptotics results (Theorems 5.1, 
5.2, 5.3, 5.4 in Section 5) on the size of the kernel, the probability of winning on trees for the first 
player... 

2. Definitions 

We give below more precise definitions, readers familiar with digraphs can skip them. 

Let D = (V, E) be a digraph. For each v G V, let = {w G V/ {v, w) G E} and v~ = {w G 
V/{w,v) G E}, \v^\ is the out degree of v and \v~\ is the in degree of v. 

A vertex with an in degree of is called a source (since one can only leave it) and a vertex 
with an out degree of is called a sink (since one cannot leave it). Let U CV, J7+ = U„g!7W+ and 
U~ = U^g[/w~, we denote by D{U) the subgraph induced by the vertices of U. 

There is a path from vertex v to w means that there exists a sequence (vi, . . . ,Vk) such that 
vi = V, Vk = w and Vi G v^_^_i U v~^i, for i = 1 ... A; — 1. There is a directed path from vertex v 
to w means that there exists a sequence {vi, . . . ,Vk) such that Vi = v, Vk = w and Vi G t^it]_i, for 
i= l...k-l. 

A cycle is a path {vi, . . . ,Vk) such that vi = Vk- A circuit is a directed path {vi,. . . ,Vk) such 
that vi = Vk- 

If D contains a directed path from vertex v tow then v is an ancestor of w and wis& descendant 
of V. If this directed path is of length 1, then the ancestor t; of w is also called a parent of w, and 

V a child of w. 

D is strongly connected if for each pair of vertices, each one is an ancestor of the other. D{U) 
is a strongly connected component of D if it is a maximal strongly connected subgraph of D. 

U is an independent set when U O — and a dominating set when OU ^ ^ for any 

V GV\U. U is & kernel if it is an independent dominating set. 

£) is a DAG if it is a directed digraph without circuit (the terminology "directed acyclic graph" 
being popular, we use the acronym DAG although it should stands for "directed acircuit graph", 
according to the above definitions of cycles and circuits). 
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3. How to find the kernel of a digraph 

Consider digraphs satisfying the following rules: 

• each vertex is colored either in red or in green, 

• each green vertex has at least a red child, 

• no red vertex has a red child. 

It is immediate to see that a digraph satisfying such coloring constraints possesses a kernel, which 
is exactly the set of its red vertices. It is now easy to see, e.g., that the circuit of length 3 has no 
kernel, that the circuit of length 4 has 2 kernels, that any DAG has exactly one kernel. For this last 
point, assume that D is a DAG (directed acircuit graph). Algorithm 1 (below) returns its unique 
kernel. It begins to color the sinks in red and then goes up toward sources, as it is deterministic 
and as it colors at least a new vertex at each iteration, this proves that each DAG has a single 
kernel. Such an algorithm was already considered by Zermelo while studying chessgame. 



Algorithm 1 The kernel of a DAG 

Input: a DAG D — (V, i?), Noncolorcd— V (i.e. no vertex is colored for yet) 
Output; the DAG, with all its vertices colored, the red vertices being its kernel 
while it remains some non colored vertices (Noncolored ^ 0) do 
for all V ^ Noncolored do 

if is a sink or if all the children of v are green then 
color V in red 

color all the parents of v in green 
remove the colored vertices from Noncolored 
end if 
end for 
end while 



For sure, it is possible to improve this algorithm by using the poset structure of a DAG, and 
thus replacing the "for all v G Noncolored" line by something like "for all v E Tocolornow" where 
Tocolornow is a set of candidates much smaller than Noncolored. 

More generally, in order to color a digraph which is not a DAG, simply split it in p components 
which are DAGs. Then, apply the above algorithm on each of these DAGs (excepted the cut 
points that you arbitrarily fix to be red or green) . It finally remains to check the global coherence 
of these colorings. As one has p cutting points (which can also be seen as p branching points in a 
backtracking version of this algorithm), this leads to at most 2^ kernels. This also suggests why 
this problem is NP: for large (dense) graph, one should need to cut at least p ^ n points, which 
leads to a 2" complexity (lower bound). 




Figure 1. The first digraph is a well-colored DAG (it has several cycles, but no 
circuit). The second digraph is a well-colored digraph (it is not a DAG, as it 
contains one circuit). The third digraph is a DAG, but is not well colored (the 
top green vertex misses a red child). [For people who are reading a black & white 
version of this article, red vertices are fulfilled and green vertices are empty circles.] 
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4. Generating functions of well-colored unicircuit digraphs 

There exists in the literature some noteworthy results on digraphs using generating functions 
(related e.g. to EGF of acychc digraphs [18, 28], Cayley graphs [26], (0,1) matrices [25], Erdos- 
Renyi random digraph model [24]), but as fas as we know we give here the first example of 
application to the kernel problem. 

The coloring constraints mentioned in Section 3 are "local" : they are defined only in function 
of each vertex and its neighbors. One nice consequence of this "local" viewpoint of kernels is that it 
opens up a whole range of possibilities for a kind of context-free grammar approach. Indeed if one 
considers rooted labeled directed trees that are well-colored {i.e. which possesses a kernel), one can 
describe/enumerate them with the help of the five following families of combinatorial structures 
(all of them being rooted labeled directed trees): 

• T: all the trees with the coloring constraint 

• TJ: well-colored trees with a red root (with an additional out-edge) 

• : well-colored trees with a red root (with an additional in-edge) 

• TJ : well-colored trees with a green root (with an additional out-edge) 

• : well-colored trees with a green root (with an additional in-edge) 

• TJ^: well-colored trees with a green root (with an additional out-edge which has to be 
attached to a red vertex) 

Those families are related by the following rules: 

Tj = X Set>i (Tj) x Set{T} U Tg^ U TJ) 
^ Tg^ =g-^ X Sct>i (Tj) X Sct{Tj: U T^^ U TJ) 
' Tj = rT X SctiTg^UTl) 
=r^ X Set(Tj UTjJ 
Jl = gT X Set(T; U U U Tj) 

The Set operator reflects the fact that one considers non planar trees, i.e. the relative order of the 
subtrees attached to a given vertex does not matter. The notation Set>i means one considers non 
empty set only. 



,1 A ,1 

AO o 
A A 

2 3 1 3 2 1 

Figure 2. A tree e Tj of size 3 and all its possible labeUings. Tj stands for 
directed trees with a green root with an additional in-edge on this root. 



As we are dealing with labeled objects (we refer to Figure 2 for the different labellings of a 
rooted directed tree), it is more convenient to use exponential generating functions, the above 
rules are then translated (see e.g. [20, 14] for a general presentation of this theory of "graphical 
enumeration/symbolic combinatorics" ) into the following set of functional equations (where z 
marks the vertices): 

'T(z) = Tj(z)+T,T(z), 

< Tj(z) = Tj{z) = z(exp(r;(z)) - 1) exp(r,i(z) + Tjiz) + Tj{z)) , 
T,I(z) = T,l(z) = zeMTjiz) + T^z)) . 

Note that T^ = T as one has the trivial bijection "Tj^ trees with a root without red child" 
= "TJ trees" and "Tj^ trees with a root with at least a red child" = "Tj trees". Define now 
Tg{z) := Tj{z) and Tr{z) :— TJ{z), the above system simplifies to: 
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'T{z) = T,{z)+Tr{z)=Tl{z), 

Tg{z) = zexp(2r(z)) - zexp(T(z) + Tg(z)) , 

Tr{z) = zcMTg{z)+T{z)) = T(z)cxp(-T,(z)) . 

This system has a unique solution, as the relations can be considered as fixed point equations 
for power series. Their Taylor expansions are: 



T{z) = z + 4^ 



2! 



36 



3! 



512 



4! 



10000 



5! 



,6 ^7 

248832— + 7529536— + O(z^) , 
6! 7! 



2 3 4 5 6 T 

Tgiz) = 2^ + 15^ + 232 + 4535|- + 114276|- + 3478083||- + 0{z^) , 

2 3 4 5 6 7 

Tr{z) = z + 2^+ 2l|j- + 280^ + 5465|j- + 134556^ + 4051453-^ + 0{z^) . 

For sure, the i-th coefficients of these series are divisible by i, as we are dealing with rooted 
object. Here are the 3 generating functions of the corresponding unrooted trees: 



T«"''-(^) = z + 2 



2! 



12 



3! 



~i ~5 _6 ^7 

128— + 2000— + 41472— + 1075648— + O(z^) , 
4! 5! 6! 7! ^ ^ 



rpunr. 

g 



58- 



z 



907- 



z 



z 



6 



19046- 



496869- 



z 



0{z'), 



z 

Tr^-{z) = ^ + 1- + 7|- + 70^ + 1093^- + 22426|^ + 578779|j- + 0{z^) . 

Of course, trees are DAG and therefore have a unique kernel. This implies that T{z) is exactly 
the exponential generating function of directed rooted trees, i.e. 

T{z) ^ C(2z)/2 and r„ = (2n)"-i 

where C{z) is the Cayley function (see Figure 3 and references [5, 10]), defined by 

n! 



C(z)=zexp(C(0)) = ^n"-i-. 




Figure 3. The Cayley tree function C{z) goes from — oo for ^ ^ — oo to 1 at 
z = -. It satisfies C{z) = 2;exp(C(2)). 



Solving the set of equations for T, Tg and finally leads to 
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Theorem 4.1 (Enumeration of well-colored trees). 
By ditrees, we mean well-colored rooted labeled directed trees. By well-colored, we mean each green 
vertex has at least a red child, each red vertex has no red child. 
The exponential generating function of ditrees is given by T{z) — C{2z)/2, 
the EGF of ditrees with a red root is given by Tr{z) = — C(— C(2z)/2), 
the EGF of ditrees with a green root is given by Tg{z) = C(2z)/2 + C(— C(2z)/2), 
where C{z) is the Cayley tree function C{z) = zex.p{C{z)). 

The EGF for the unrooted equivalent objects can be expressed in terms of the rooted ones: 
yunr. =T -T^ , T'g"'''- = T™''- - r;™''- , and T™''- = 2T - 2TTr + Tr- 2T/Tr + T^ /2 . 

Proof. The formulae for T,Tr and Tg can be checked using the definition of C{z) in the 
fix-point equations in the simplified system above. The fact that the GF for unrooted trees can be 
expressed in terms of the GF of rooted ones can be proven by integration of the Cayley function, 
or by a combinatorial splitting argument on trees. □ 

We can go on and enumerate the different possibilities of circuits for a well-colored digraph. 
They can be described as 

Gyc(g) U Cyc(r -^{g ^}+) 
This reflects the fact that either a circuit is made up of green vertices only, or it contains some 
red vertices, but they have to be followed by at least a green vertex. NB: Whether one counts or 
not the cycles of length 1 (i.e. a single red or green vertex) will only modify the first term of the 
generating function. Symbolic combinatorics [14] translates the above cycle decompositions in the 
following function: 

In I + In ' ^ ^ 



l-gj ' 

y / \ l-g / 

where r / g mark the number of red/green vertices. This leads to the following Theorem: 

Theorem 4.2 (Enumeration of possible well-colored circuits). 
The exponential generating function of possible well-colored circuits is given by 

L{z) = - ln(l ~z-z^) = z + ?,^+^^+ 42^7 + 264^- + 2im^ + 2088o|^ + O(z^) . 

2! d! 4! 5! o! 7! 

Its coefficients satisfy i„ ~ (n — 1)! (0" + (1 — 0)"), where i„ are known as the n-th Lucas number 

(usually defined by the recurrence Ln+2 ~ L^+i + Ln,Li = 1,L2 = 3) and where (J) — {1 + \/5)/2 

is the golden ratio. 

Note that a reverse engineering lecture of this generating function leads to the simpler de- 
composition Cyc((7 U rg), which also explains the recurrence! Now, the following decomposition of 
possible cycles is trivially related to the decomposition of possible circuits: 

Cyc(r x{^gU ^g}+ x { ^ U ^}) U Cyc(5 ^ U <? ^) 

leads to the EGF — ln(l — 2z — 4z^) whose coefficients are, with no surprise, 2"L„. 




_To- 



Tct 



Figure 4. Unicircuit digraphs consist in a circuit with attached trees on it. The 
left picture above is a unicircuit digraph, to the right, we give its "canonical 
decomposition" as a circuit of atoms which are trees. Any well-colored unicircuit 
digraph has such a "canonical decomposition" . 
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Using the decomposition given in Figure 4, one obtains the generating function for unicircuits; 

Theorem 4.3 (Enumeration of unicircuit well-colored digraphs). 
The EGF of unicircuit well-colored digraphs is 

U{z) = T""''- - Tg + In ' ^ 



1 - (Tg + TgTr) 



C{2zf _ c^_C(^^ _ In - ^ - + C{- ^) ^ 



4 '2' V 2 ^2' ^ 2 

= z + 4|j- + 30|| + 452^1- + 8840^ + 224832^ + 6909784||- + 0{z^) , 

where C{z) is the Cayley tree function C(z) = 2:exp(C(2;)). 

Now, consider the larger class of unicycle digraphs (digraphs which have or 1 cycle). Recall 
that a circuit is a cycle, but a cycle is not necessarily a circuit. In order to get a "canonical 
decomposition" for unicyle digraphs (similar to the one given in Fig. 4 for unicircuit digraphs), one 
considers 3 cases: 

• Either the graph has no cycle, those graphs are counted by T""''- . 

• Either it is a cycle with only Tg trees branched on it (i.e. no red nodes in the cycle), those 
graphs are counted by (In (^ i_\t ) ^ — 4Tg^/2)/2 + where 2Tg corresponds 
to Tg X {— > U ^}, one removes cycles of length 1 and 2 from the logarithm (this explains 
the —2Tg — ATg'^/2 term) and one divides the whole formula by 2 because one has to 
take into account the fact the cycle can be read clockwise or not, and one adds the only 
legal cycle of length 2. 

• Either the graph contains a cycle with some red nodes and then one considers the following 
possible "bricks": 



(T 4- 










Tr^ 






(but not 


a cycle of length 2, because multiarcs are not allowed) 


Tr- 


- {TA- 


^ U 4 


-}rTg^ 


(but not a cycle of length 2) 


' T - 


- (TA- 


+ U 4 


-}rTg^- 




T <- 


-T,A- 


^ U 4 


-}{Tg{^ 


u-})*r,4- 


Tr- 




^ U 4 


-}{Tg{^ 


u -})* Tg^ ^ 



Theorem 4.4 (Enumeration of unicycle well-colored digraphs). 

The EGF of unicycle well-colored digraphs is 

V{z) = T^^^- + lln{^^^~Tg-Tgy2-TrTg/2-TrT/2 




1-2T„ 



ln(l + T^) - ^ ln(l - 2T) 

= z + A^+ 36^ + 692^ + 15920^ + 458622^ + 15559264^ + Oiz^) . 
2! 3! 4! 5! 6! 7! 

where T, Tg, Tr, and T""""- are given in Theorem 4-1- 

Note that in the two theorems above, any given non-colored graph is counted with multiplicity 
0, 1 or 2 (if there are 0, 1 or 2 ways to color it). We explained in Section 3 that a multiplicity 
larger than 2 was not possible for unicycle digraphs. We enumerate in the following proposition 
those with exactly 2 possible colorations. 

Proposition 4.5 (Enumeration of unicycle digraphs with two kernels). 
The EGF of unicycle digraphs with 2 kernels is 

D{z) = -ln^/l + C{-C{2z)/2)\ 

where C{z) is the Cayley tree function C{z) = zexp{C{z)). 
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Remark: From the definition of cycle/circuit, D(z) is also ttie EGF of unicircuit digraphs with 
2 kernels. 

Proof. Let T) be the set of unicycle digraphs with 2 kernels. First, it is easy to see that 
Cyc(T^) C V (with a sligth abuse of notation, as wc first only consider the shape, not the coloration 
of the Tj. trees): simply color the nodes in the cycle alternatively in green and red, and switch the 
colors of a part of attached trees, if needs be. 

Wc now prove the next step V C Cyc(T^): Take a unicycle graph in I?, it means at least one 
of its vertex can be colored both green and red. Such a vertex v can be taken, without loss of 
generality, in the circuit (from the above remark, the cycle is in fact a circuit). [If it were not the 
case, all bi-colorabled vertices would be in the tree components, but then one could split our graph 
to get DAGs which are known to be uniquely colorable]. But when v is red, it implies it has only 
Tg trees attached to it, which means than when it gets green, the next node in the circuit has be 
red (and was previously green!). This implies alternation red/green (and even length for parity 
reasons) for all the nodes in the circuit. 

This leads to a canonical decomposition 

If one divides by 2 for the (anti)clockwise readings, this leads to the Theorem. □ 

Most of these results (and also the computations of Section 5 hereafter) were checked with 
the computer algebra system Maple. A worksheet corresponding to this article is available at 
http://algo.inria.fr/banderier/Paper/kernels.mws (or kernels .html), it uses the Algolib 
librairy, downloadable at http://algo.inria.fr/libraries/). 

5. Asymptotics 

In this section, we give asymptotic results for n +oo. 

Theorem 5.1 (Proportion of trees with a green/red root). 
AsymptofAcaMy w 47.95% of the trees have a green root, where the constant A « 0.351733 is 
defined as the unique real root o/2A = exp(— A). 

A more pleasant way to formulate this Theorem consists in considering Nim-type games (first 
player who cannot move loses) on directed, trees where the tree and the starting position are chosen 
uniformly at random. The strategies of the two players being optimal, the first player has then a 
probability of 4'^. 95% (asymptotically) to win the game. (Recall that if the starting position can he 
chosen by the first player, then he will always win.) 

Proof. The key step of this result and the following ones are the following expansions (derived 
from the expansion of the Cayley function) for T, Tr and Tg-. 

T{z) ~ ^ - A=Vl - 2ez + 0(1 - 2ez) 
6 v2 

Tr{z) - A - -^Vl - 2ez + (9(1 - 2ez) 
1 + A 

T,(z) ~ 1 - A - ±i-^VT^^z + 0(1 - 2ez) , 

where the constant A is defined as A := Tr{^) ~ 0.351733. 

By Pringsheim theorem [14], as Tr{z) has nonnegative coefficients, then Tr{z) has a positive 
singularity. As coefficients of Tr are smaller than coefficients of T, its radius of convergence belongs 
to [0, l/(2e)]. Now, —C{2z)/2 is negative on this interval, and thus C{—C{2z)/2) is analytic there, 
and l/(2e) is therefore its only possible dominating singularity. The radius of Tg follows from 

T = Tr + Tg. The theorem follows by considering 'ffj'^^g = TTX " ^ + ° 

Theorem 5.2 (Proportion of red vertices in possible circuits). 
Asymptotically \ — —7= ~ 27.63% of the vertices of a possible circuits are red. 
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Proof. One has to considcrer the fohowing bivariate generating function (exponential in 
z, ordinary in u): In ^ i-(zXuz'^) )- "^^^ wanted proportion is then given by ^^jJqf^Tij^i where 
[z"']duF{z, 1) means the n-th coefBcient of "the derivative with respect to u of F{z,u), then eval- 
uated at u = 1" . □ 

Then, one can wonder if the asymptotic density of well-colored unicircuit graphs is more than 
50% or even if it is 100%? The following theorem gives the answer: 

Theorem 5.3 (Proportion of well-colored unicircuit digraphs). 
The proportion of well- colored graphs amongst unicircuit digraphs is asymptotically: 

3A3 + A2 - A - 1 



(1 + A)2(A-1)) 
where A is the constant defined in Theorem 5.1. 



92.65% 



Proof. Relies on a singularity analysis of the generating function of Theorem 4.3, with the 
expansions given in Theorem 5.1. Note that some unicircuit digraphs can have 2 kernels, so one 
has to perform the following asymptotic expansions: 

™4^. 92.65-^ + 0(J,), 

[z-^]F{z) n ^n^' 

where F{z) = T^^''{z) + ^'a{ i_T(z) ) ~ ^(•^) the EGF of (non-colored) unicircuit digraphs. □ 

For sure, it one considers now the asymptotic density of well-colored unicircuit graphs, the 
proportion should be larger, as one only adds DAGs (which are all well-colorable). The following 
theorem gives the noteworthy result that unicircuit graphs are in fact almost surely well-colored: 

Theorem 5.4 (Proportion of well-colored unicycle digraphs). 
There is asymptotically a proportion of 1 — ^ « 1 — of well-colored graphs 

amongst unicycle digraphs of size n, where A is the constant defined in Theorem 5.1. 

Proof. Relies on a singularity analysis of the generating function of Theorem 4.4, with the 
expansions given in Theorem 5.1. Note that some unicycle digraphs can have 2 kernels, so one has 
to consider 

[z'']V{z) - D{z) 

where G{z) = r""'-(z) + 5 111(33^^) - T{z) - T{zf/2 is the EGF of (non-colored) unicycle 
digraphs (one substracts T^/2 because amongst the 4 graphs with a cycle of length 2 created by 
the ln( i_2T{z) ) P^rt, 3 are not legal: 1 was already counted because of symmetries, and the other 
2 have in fact a multiple arc, whereas it is forbidden for our digraphs). □ 

Finally, if one considers graphs with at most k cycles, it means one has more cutting points, 
which relaxes the constraints for well-colarability (=existence of kernel). According to the above 
results, this implies an asymptotic density of one. This gives as a corolary of our results, that all 
these families have almost surely a kernel. A kind of "limit case" is dense graphs, for which some 
results already mentionned by Fernandez de la Vega [13] and Tomescu [30] give that they have 
indeed almost surely a kernel. 
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6. Conclusion 

It is quite pleasant that our generating function approach allows to get new results on the kernel 
problem, giving e.g. the proportion of graphs satisfying a given property, and new informations on 
Nim-type games for some families of graphs. 

As a first extension of our work, it is possible to apply classical techniques from analytic 
combinatorics [14] in order to get informations on standard deviation, higher moments, and limit 
laws for statistics studied in Section 5. 

Another extension is to get closed form formulas for bicircuit /bicycles digraphs, (the generating 
function involves the derivative of the product of two logs and the asymptotics are performed like 
in our Section 5). It is still possible (for sure with the help of a computer algebra system) to do 
it for 3 or 4 cycles but the "canonical decompositions" and the computations get cumbersome. In 
order to go on our analysis far beyond low-cyclic graphs, one needs an equation similar to the one 
given by E.M. Wright [32, 33] for graphs. Let Wt be the family of well-colored digraphs with 
£ edges more than vertices, (£ > —1). It is possible to get an equation "a la Wright" for We by 
pointing any edge (except edges linking a green vertex to a red one) in a well-colored digraph. It is 
however not clear for yet if and how such equations can be simplified in order to get a recurrence 
as "simple/nice" to the one that Wright got for graphs, thus opening the door to asymptotics and 
threshold analysis beyond the unicyclic case. 
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